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ABSTRACT 


There are many fluid flows where the onset of transition can be caused by dif- 
ferent instability mechanisms which compete among themselves. Here we consider 
the interaction of two types of instability mode (at an asymptotically large Reynolds 
number) which can occur in the flow above a rotating disc. In particular, we exam- 
ine the interaction between lower-branch Tollmien-Schlichting (TS) waves and the 
upper-branch, stationary, inviscid crossflow vortex whose asymptotic structure has 
been described by Hall (1986). This problem is studied in the context of investi- 
gating the effect of the vortex on the stability characteristics of a small TS wave. 
Essentially, it is found that the primary effect is felt through the modification to 
the mean flow induced by the presence of the vortex. Initially, the TS wave is taken 
to be linear in character and we show (for the cases of both a linear and a nonlinear 
stationary vortex) that the vortex can exhibit both stabilizing and destabilizing 
effects on the TS wave and the nature of this influence is wholly dependent upon 
the orientation of this latter instability. Further, we examine the problem with a 
larger TS wave, whose size is chosen so as to ensure that this mode is nonlinear in 
its own right. An amplitude equation for the evolution of the TS wave is derived 
which admits solutions corresponding to finite amplitude, stable, traveling waves. 


* This research was supported by the National Aeronautics and Space Adminis- 
tration under NASA Contract No. NAS 1-18605 while the author was in residence 
at ICASE, NASA Langley Research Center, Hampton, VA 23665. 



1. Introduction 


Many studies concerned with the instability of three-dimensional boundary 
layers have been motivated by a desire to understand the phenomenon of transition 
to turbulence in fluid flows. Here we are concerned with a self-consistent asymptotic 
description of the interaction of stationary cross-flow vortices and lower-branch 
Tollmien-S chlichting waves (hereafter referred to as TS waves) in the boundary 
layer of the flow induced by a rotating disc. This particular flow is susceptible to an 
instability similar to that which occurs in the boundary layer of flows over a swept 
wing; a situation which has practical relevance to the development of Laminar Flow 
Control wings. Further, there is an exact solution of the Navier-Stokes equations 
which describes the rotating disc flow. This makes the study of this rotating disc 
problem particularly suitable for a theoretical analysis of the interaction of the 
vortices and the TS waves. 

We will concentrate on the description of the interaction which occurs in flows 
at an asymptotically large Reynolds number (which is based upon the angular 
velocity of the disc, a typical lengthscale of the problem and the kinematic viscosity 
of the flow). The lower-branch TS waves are then described by a classical interactive 
triple deck structure, the key elements of which are reviewed by Messiter (1979), 
Stewartson (1981) and Smith (1982). The TS waves, which are travelling modes, 
have a wavelength much greater than the boundary layer thickness and a small 
wavespeed. The structure of these disturbances is outlined in Section 3. 

The crossflow vortex instability structure can occur only in three-dimensional 
boundary layers and was first examined both theoretically and experimentally by 
Gregory et al (1955). The stationary vortex mechanism appears when a direction 
for the disturbance is chosen such that the effective basic velocity profile contains 
an inflexion point at the same location at which it vanishes. Gregory et al used 
the china clay technique to show that for the rotating disc flow the vortex insta- 
bility takes the form of a regularly spaced pattern of equiangular spiral vortices 
which is stationary relative to the disc. Stuart (in Gregory et al), using inviscid 
theory, predicted the number of vortices which would be observed in experiment as 
approximately four times greater than that actually seen, although his calculation 
of an angle of 13° between the axes of the vortices and the radius vector on the 
disc was in excellent agreement with the experiments. The difference between the 
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experimental observations and the results of the inviscid analysis has been shown to 
be due to viscous effects. Malik (1986a) calculated the neutral curve for stationary 
disturbances and found a second mode. Earlier, Federov et al (1976) had observed 
experimentally this second mode, which, as in the inviscid case, appeared as a pat- 
tern of spiral vortices. The number of these vortices was seen to be between 14 and 
16 and these had axes inclined at angles of approximately 20° to the radius vector. 

Hall (1986) has given a linear, asymptotic account of the inviscid mode found 
by Gregory et al (1955) for large Reynolds numbers. This work has been extended 
by Gajjar (1989) to examine nonlinear effects. Hall (1986) also elucidated a triple- 
deck type structure for the second type of stationary vortex which corresponds to 
an effective velocity profile with zero shear stress at the wall. This problem has 
also been studied using a weakly nonlinear approach by MacKerrell (1987, 1988). 
Further, Bassom & Gajjar (1988) have examined the properties of a non-stationary 
version of the crossflow instability. 

Much work has been performed in relation to the important problem of the 
interaction of the interaction of TS waves with Gortler vortices (an instability as- 
sociated with flows over curved surfaces). See, for example, Nayfeh (1981), Malik 
(19866), Bennett & Hall (1988), Bennett et. al. (1988), Hall Smith (1988, 

1989a, 6), Daudpota et al (1988), Bassom &; Hall (1988) and the references therein. 
However, relatively little attention has been paid to the type of crossflow-TS in- 
teraction examined in this paper, although, in particular, we refer to the work of 
Reed (1984, 1985) who investigated this interaction and found that the crossflow 
vortices lead to a ‘double exponential’ growth in the TS waves. A linear analysis of 
the problem was used and the vortex was permitted to grow exponentially although 
no account was made for nonlinear effects. The vortex was allowed to force a TS 
wave of appropriate wavelength such that the growth rate of the latter mode was 
proportional to the amplitude of the vortex. The amplitude of the TS wave then 
grows like the exponential of an exponential. In the aforementioned papers Reed 
also examines the question of interaction between crossflow vortices of particular 
wavelengths. However, in both of these problems, the approach adopted ignores 
the crucial fact that for a completely rational description of the importance of the 
interaction process, each of the instability modes involved should be neutrally sta- 
ble at leading order in their own rights. Otherwise, the vortex growth due to any 
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interaction process is no larger than the growth experienced in its’ absence and the 
importance of the role of the interaction is impossible to assess. For this reason, 
together with the neglect of nonlinear terms in her subsequent analysis, it must 
be concluded that the relevance of Reed’s work to practical situations is at best 

doubtful. 

The primary motivation for the present investigation stems from the question of 
how the presence of an inviscid stationary vortex affects the stability characteristics 
of the TS wave. We emphasize that at this stage we do not concern ourselves 
with calculations for the growth rates of the respective instabilities. Firstly, we 
consider the problem of a small crossflow vortex which has an asymptotic structure 
as described in Hall (1986). The vortex has a critical layer structure in the vicinity 
of the inflexion point of the basic flow and this critical layer is linear in character. 
The vortex then has a very small TS wave superimposed upon it and the effect 
of the vortex on the neutral stability properties of the TS wave is obtained. We 
find that the vortex can stabilise or destabilise the TS wave depending upon the 
orientation of the latter mode. The amplitude of the TS wave is then allowed to 
increase to the point at which it becomes nonlinear in its own right and amplitude 
equations for the evolution of the TS wave are obtained. 

In addition, the problem of the interaction between a stronger, nonlinear vortex 
and a small TS wave is considered. It is concluded that for both the linear and the 
nonlinear vortex problems the effect of the vortex on the stability of the TS wave 
is felt primarily through the correction to the effective mean flow induced by the 
presence of the vortex. 

The procedure in the remainder of the paper is follows. In Section 2 we formu- 
late the problem and indicate the basic disturbance structure for the linear crossflow 
vortex. We impose a small TS disturbance on this configuration in Section 3 and 
we examine the neutral stability characteristics of this wave in the presence of the 
vortex. The TS wave is increased in size to become nonlinear in Section 4 and the 
problem with a larger, now nonlinear, vortex is considered in Section 5. Finally, we 
draw some conclusions in Section 6. 


3 



2. Formulation of the problem and the linear vortex structure 

We consider the case in which the disc rotates about the *-axis with angular 
velocity Cl. Relative to cylindrical polar axes ( r,6,z ) which rotate with the disc 
and in which r and z have been made dimensionless with respect to some reference 
lengthscale l, the continuity and Navier-Stokes equations for an incompressible fluid 
in the region z > 0 are 

V.u = 0, (2.1a) 

— + (u. V) u + 2 (n A u) + n A (n A r) = - i Vp + u V 2 u. (2.16) 

Here u is the velocity vector, r is the coordinate vector, p is the fluid density, p is 
the fluid pressure and u the kinematic viscosity. The Reynolds number R for the 
flow is given by R = Cll 2 /v and is taken to be large throughout the following. 

Anticipating the structures which govern the interactions described in the fol- 
lowing sections it is found convenient to define the small parameter e = R~t*. As 
the axes rotate with the disc, the basic steady flow is given by the Von-Kdrm^n 
solution 


u = u B = in (ru(£),rv(£),e 12 tD(£)) , p = p/ 2 fi 2 e 24 p(f), (2.2) 

where z = e 12 f and u,v,u> satisfy the equations 


— 2 r - >2 1 n 

« — (1 + V) + u w — u , 


2u (1 + v) + v w = v 


(2.3a, 6) 


w + 2u = 0, 


_> i _« 

p + to to — to = 0. 


(2.3c, d) 


Here primes denote differentiation with respect to £ and the appropriate boundary 
conditions are 


u=v = to = 0 on £ = 0, 
u — ► 0, v — ► —1 as £ — ► oo. 


(2.3e) 


We shall be concerned with perturbations to the basic flow and it is convenient 
to now derive the equations governing these disturbances. If the steady solution is 
perturbed by writing 


u = u B + ICl (U, V,W), p = pl 2 Cl 2 (p + P) , (2.4) 
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where U,V,W,P are small three-dimensional disturbances and if expressions (2.4) 
are substituted into the governing equations (2.1) for the flow in the rotating frame, 
we then obtain the following perturbation equations. 


U dU ldV dW_ 
r ^ dr r dO dz 

_ _ 2 

— _j_ r u— + v-^r + e 12 w-^~) U + uU + rW — - 2(1 + v)V — 
dt dr dd dz) dz r 

(4 + 7Te + w h) u = + H m " ? " ?&) ’ 

l +r 4 + 4 +£l2 4) v+EV+rW, S +2(1+ ^ + ^ 

( d Vd d\ 1 dP 1 V 2 dU\ 

+ { u l* + 7Te + w ^) V= -TM + R\ LlV) -^ + ^M)' 

(s + s + a +( , 2 a\ 

\dt dr d9 dz J 


( 


+ 


( 


(■ 






t0 „dw 
W + e 12 ^ — 
dz 


Here, 


a 2 i a 2 i a_ a 2 

L= dr 2 + T 2 d9 2 + r dr^ dz 2 ' 


(2.5a) 


(2.56) 


(2.5c) 


(2.5d) 


2.1 The linear vortex structure 


Since we will be concerned with the effect of the vortex on a small TS wave 
which, in the first instance, will be taken to be of such size so as to have an in- 
finitesimal influence on the vortex, we first need to consider the structure of the 
vortex when the TS wave is absent. Following the work of Gregory et al (1955), the 
inviscid vortex has wavelengths scaled on the boundary layer thickness and so we 
consider disturbance quantities with r and 6 dependence given by E, defined by 


E = exp 


a(r, 4 )* + W(t)} . 


(2.6a) 


We shall restrict our attention to examining disturbances in a neighbourhood of 
some point r = r n and expand a and /? in the forms 


a = <* 0 4- e 4 ai + . . . , P = Po + c 4 /?i + (2.66) 
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where 7o — a o + 1S the effective wavenumber of the disturbance. Thus W 10 
satisfies a Rayleigh equation and we obtain stationary modes by demanding that 
the point of inflexion of the effective basic flow coincides with the point at which this 
flow vanishes, say at £ = £. This criterion, together with ( 2 . 9 ) and its appropriate 
boundary conditions that W 10 — > 0 as £ — ► 0 and as £ — * oo yields the values 


for otQ and (3q. The solution of this eigenproblem gives 

i = 1 . 46 , = 1 . 16 , ^ ( 2 . 10 ) 

Po r n 

In this inviscid zone it is also found that the mean flow quantities U 30 , V 30 , W 30 , 
P 30 satisfy the equations 

217 30 + W zo = 0, (2.11a) 

^30 ~ WU 30 ~ 2fit/ 3 o — u W 3 0 + 2(1 + = 0, (2.116) 

^30 — — 2uV 3 o — v W 3 0 — 2(1 + v)U 3 0 = 0, (2.11c) 

P 30 = -2 (l^ol 2 ) • (2. lid) 


To determine these correction terms we need to consider explicitly both the 
wall layer and the critical layer structure (at £ = £). In the wall layer where we 
define the O(l) co-ordinate Z by £ = e 4 Z , the total velocity and pressure fields 
assume the forms 


u — r n (u ) Z + . . . + 6 [(I/*, . . ,) 2 ? + c.c.j 

+S 2 [(C/Jo + ...)E 2 + c.c. + r n e~ 8 (t4 + •••)]+..., 
v = ^n(v )s=q£ 4 Z + . . . -f 6 J(4) + . . .)E + c.c.j 
+6 2 [(4 + . . .)E 2 + c.c. + r n c- 8 (4 + ...)]+..., 


w — 2^ W )$=o e2 °^ 2 + • • • 


+ 6 


e*(wl + ...)E + cx. 


+ e4 ( W 20 + • • -) e2 + c.c. + e 8 (l 4 + •••)] + ••• , 

P = e 24 (p)^=o + ■•• + £ e 4 (p} 0 + . . .)E + c.c.j 

e 4 {P ^ 0 + ...)E 2 + c.c. + (pj 0 + ...)]+.... 


+ 6 X 


( 2 . 12 a) 


(2.126) 


( 2 . 12 c) 


(2.12d) 
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In general, the wavenumbers will be complex quantities but we shall only con- 
sider neutral disturbances and so take ao,a u ...,l3o,Pi--- to be real. The dis- 
turbance structure in the 2-direction is as described in Hall (1986), Gajjar (1989). 
There is an inviscid layer of 0(e 12 ), (the same depth as the boundary layer), and 
to satisfy no-slip conditions at the wall a viscous layer of thickness 0(e 16 ) must be 
present. In contrast to the earlier work we need to explicitly consider the mean flow 
correction induced by the vortex because this correction plays a crucial role in the 
interaction problems to be described in subsequent sections. We find that in the 
inviscid zone the velocity and pressure fields assume the forms 


u = r n u(£)+6 [(U 1Q + t A U\i + . . .) E + c.c.] 

+ S 2 [(17*0 + ...)E 2 + c.c. + r n e~ 12 (U 30 + ...)]+ 

v = r n i>(£)+<5 [(Vko + + • • •) E + c.c.] 

+ <5 2 [(V20 + . • •) E 2 + c.c. + r n e -12 (V30 + •••)]+ 0(6 3 ), 

w = e 12 u)(^)+5 [(W 10 + e 4 W n + . . .) E + c.c.] 

+ 8 2 [{W 20 + ...) E 2 + c.c. + (W 30 + ...)]+ 0{6 3 ), 

p — e 24 p(£)+^ [ (-Pio + ^P\ \ + ..•)■£ + c - c - ] 

+ 6 2 [(P 20 + • • ■) E 2 + c.c. + (P 30 + •••)]+ <>(6 3 ). 

Here 6 *C 1 is the infinitesimal vortex amplitude, the disturbance quantities 
U 10 ,U 1U . . . , U 2 o, U 3 o etc. are functions of £ and r n and c.c. denotes complex con- 
jugate. We substitute (2.7) into (2.5) and find that the fundamental vortex terms 

satisfy 


(2.7a) 

(2.76) 

(2.7c) 

(2.7<f) 


i ^aot/io + + ^10 — Oj 

(2.8a) 

iU OB Uio + r n u W 10 = -ia 0 Pio, 

(2.86) 

t , fio 

iUoB^lO + TnV W10 = — * P\ 0 i 

r n 

(2.8c) 

/ 

iUosWio = —P io> 

(2.8d) 


where U 0B = r n a 0 u + (3 0 v is the effective basic velocity profile. Eliminating Uio,V i0 
and Pio from these equations yields 


UoB ( w ','o - lo^io) = 'C B W U> , 


( 2 . 8 ) 
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Substituting these expansions into (2.5) shows that wf Q ,pf 0 satisfy the 

equations given in Hall (1986) for the wall layer quantities. The solution of these 
equations ultimately leads to the determination of the values of a\ and (3\ in 
(2.6). However, this is not of primary interest here. We remark that the quantities 
^30 > V&, ^30 are proportional to Z and respectively, and hence the boundary 
conditions at £ = 0 for the mean-flow disturbance quantities f/ 30 , V30, W 30 (given 
by (2.11) above) are that 

^30 = ^30 = W30 = 0 on £ = 0. (2.13) 

Clearly, we also require that the mean flow quantities decay as we leave the boundary 
layer so 

Uzo,V 30 — .0 as £ — >oo. (2.14) 

However, we also need to determine the effect of the critical layer at £ = £ 
on these mean flow terms. From (2.8) we see that as £ — > (, U 0B = 0(£ - £) 
and singularities exist in the inviscid zone solutions U 10 and Vx 0 . To smooth these 
singularities we invoke viscous effects by considering a zone of thickness O(e 10 ) 
surrounding £ = £. In this region we suppose that £ = £ + e 4 z where z is 0(1) and 
then the velocity and pressure fields are 

« = ^r n U| + e 4 r n u'gz + ^e 8 r n u'^z 2 + . . . j + 

^[( e 4 ^10 + £^11 + e 4 [7i2 + . ..)£? + c.c] 

+ <5 2 [r n e 12 ( Uzo + £ 4 f?3i + 

v = + e 4 r n v'^z + ^e 8 r n v'^z 2 + . . + 

6 [(e~ 4 Vio + V n + £ 4 F 12 + ...)£ + c.c.] 

+ [ r n c 12 (^30 + t 4 ^31 

w = c 12 (w£ + e 4 tU|Z + . ..) + 6 [(W 10 + e 4 Wn + e 8 W X2 + . . .) E + c.c. 

+ 6 2 [(W 30 + e 4 Wz l + . ..) + ...], 



(2.15a) 


(2.156) 


P — e 24 p(0 + . . . + 6 [(P10 + e 4 Pn + e 8 P 12 + . . .) E + c.c.] 
+ 6 2 [(-P30 + Pz\ + ...) + ...]. 


(2.15d) 
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The solutions for the fundamental terms are well known. Substituting these 
expansions into (2.5) and recalling the definitions of Uio,V lo ,W lo and Pio we find 

that t J5 q i 

Wio = VM£), -Pio = Pio(e) = —r-^ioCC), 

/o 

where we have defined Bjk = r n ajU^ k ^ (£) + (3 jV^ H 0 • 

iFrom continuity we find that V 10 = -^^io and equation (2.56) implies 

that 


ip 

dz 2 


— i(Bo\z + B\o)Uxo 


-^■(/?0«£ - r n O£OU|)Wlo(^)- 
r nl 0 


(2.16) 


To match with the solutions in the inviscid zone requires that U io - 0(1/ z) as 
|2| ► oo and the solution of (2.16) with this property is 

■tijjK'.od) j exp [ir (z + §^)] exp ( 5^) dr, 

(2.17) 


fio 


UlO = — - 

r n lo B 01 * 


r n a 0 v>) 


where we remark that Boi < 0* 

Turning to the equations for the mean flow terms in (2.15) we find that W zo is 
constant across the critical layer, but that Uzo satisfies 


dU 30 

dz 


= - + n'.off.-o) . 

r n 


or , - 

Vzo = — / R‘(0 w w; 0 )dt + A, (2-18) 

r n J — oo 

where A is some constant and an asterisk on a quantity denotes the complex conju- 
gate of that quantity. It is apparent from (2.18) that, unless Re(U l0 W/ 0 )dt = 0, 
the mean flow term C/ 30 in the inviscid zone must suffer a discontinuity across the 
critical layer. We may determine this jump by multiplying (2.17) by Wf 0 , inte- 
grating with respect to z and interchanging the order of integration. The resulting 

integral assumes a form of the type 

J{i)= j~ M{t) (ciz±y t , (2.i9) 

and the behaviour of J(i) for large |z| is well known, see Haberman (1976). After 
some algebra, taking real parts and inserting in (2.18), we retrieve the necessary 
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discontinuity across the critical layer, 


[0ao(*)]“ 


oo 

oo 


2irPo(p 0 u'; - r n a 0 V;l 


(2.20a) 


We also find from the governing equations that V 30 = -r n a 0 U 30 /Po so that 


[M*)l 


oo 
— CO 


fn&O 

00 


p3o(s)]r 


oo 

oo 


(2.206) 


Moving on to higher order terms in (2.15a - d) reveals that the term dU 3i /dz 
also suffers a non-zero jump across the critical layer which means a discontinuity 
in the derivative of the mean flow term U 3 q (in the main boundary layer) across 
£ = £• Repeating the procedure described above we find that 


= Pso]!**,, (2.20c) 

— CO 

and 

CO 

= ®|[^3o]!° 00 . (2.20 d) 

— OO 

To conclude this section, we note that the presence of the vortex has induced a 
mean flow correction term of size 0(6^e ^) across the whole of the boundary layer. 
This correction is governed by equations (2.11) with boundary conditions (2.13), 
(2.14), and jump conditions hold across £ — £, which are given by 


dV 31 

dz 


dUz i 

dz 


[£M£)]ft = 
[v»({)l|! = 


2tt/?o(/W? - r n a 0 i4) 


01 


r n°L 0 

Po 


lo r lB t 
[^30(f)]fc, 


|1Vio(£)| : 


(2.21a, 6) 


'dU, 


30 




L d(L J 


= ® { -[V30||l. 


f- 


'dV, 


30 


ii+ 


L dt 


= ^-[^3oj|! 


(2.21c, d) 


We have concentrated on the properties of the mean flow correction terms be- 
cause when we consider the interaction of the vortex with the TS wave in subsequent 
sections, it is found that the effect of the vortex on the XS wave is felt primarily 
through this correction to the mean flow. 

We solved the differential system (2.11), (2.13), (2.14), (2.21) using a fourth 
order Runge-Kutta method combined with a shooting technique. Of importance 
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is the sign of the discontinuity appearing in (2.21a) and it was found that this 
term is negative. Clearly, the scaled size of the vortex affects the magnitude of 
the discontinuities in (2.21) and in Figure (l) we illustrate the solutions for the 

corrections C/ 30,^30 and W 30 when [t/ 3 0 (€)]|! = - 1 - We find in this CaSe that 
t/' 0 (£ = 0) = -0.86 x 10- 2 and F30U = °) = “0.896. We now consider the effect 

of adding a small TS wave into the flow. 


3. Interaction of a linear vortex with a small TS wave 


The lower-branch TS waves we shall concern ourselves with here are the classi- 
cal type described by a triple deck structure. The waves have 0(e 9 ) wavenumbers 
(and hence wavelengths much larger than the 0(e 12 ) thickness of the boundary 
layer) and small, 0(e 3 ), wavespeeds. The details of the triple deck may be found 
in, for example, Smith (1979a), but we note here that the three decks are of thick- 
nesses 0(e 9 ), 0(e 12 ) and 0(e 15 ) in the direction normal to the disc. The first of 
these decks is a region of potential flow, the middle deck is primarily inviscid but 
rotational in character and the thin wall layer is viscous. 

We shall suppose in the first instance that the size of the TS wave is extremely 
small, say A(< 6) where we recall that <5 is the infinitesimal vortex size. In this case, 
at leading orders at least, the vortex structure described in the previous section is 
unaffected by the presence of the TS wave. We have already shown that the vortex 
induces a correction of size O (6 2 R?) to the mean flow throughout the boundary 
layer and we anticipate that there are essentially two mechanisms which could be 
responsible for altering the neutral stability characteristics of the TS wave. Firstly, 
we have the modification to the mean flow due to the vortex and, secondly, the 
TS wave will interact with the vortex terms to produce further modes which in 
turn interact to affect the TS and vortex instabilities. Of crucial importance for the 
ensuing study is the question as to which of these processes is the more important as 
far as the TS wave is concerned. We find after analysis that the former mechanism 
(the mean flow change induced by the vortex) has the dominating effect and this 
simplifies the following calculations. 

As the vortex induces a relative change of 0(6 2 e -12 ) in the size of the mean 
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flow we expect a similar effect on the neutral wavenumbers and wavespeed for the 
TS wave. Anticipating this result we seek TS disturbance quantities with r and 0 
dependence given by F , defined by 

F = ex P £{/ “Ke)dr + ^( e )-n(e)t|j, (3.1a) 

where, 

a = (a 0 + ...) + S 2 e~ 12 (a 00 + ...), + ...)+ 6 2 e~ 12 {0 oo + ...), 

n = e 3 {n 0 + ...)+ <s 2 e~ 9 (n 0 o + ...). 

(3.16) 

We see that the perturbed wavenumbers a 0 o,/?oo and the perturbed frequency f2 00 
measure the effect of the vortex on the neutral stability of the TS wave. Of course 
there are terms in the wavenumber and frequency expansions not explicitly detailed 
in (3.1) and larger than the a oo>/?oo a °d floo corrections, but these terms are 
independent of the vortex properties whereas our aim here is to determine this 
vortex effect. Hence we shall merely seek expressions for the values of aoo,/?oo and 
ftoo- 

The fundamental vortex terms and the TS wave interact to produce insta- 
bilities with spatial dependences of the forms EF, E~ l F etc. These waves have 
a structure similar to that for the crossflow vortex but with a more complicated 
critical layer configuration surrounding £ = £. Indeed, here these waves have two 
critical layers separated by an 0(e 1B ) distance and each of thickness 0(e 16 ). There 
is the potential for the interaction of these modes with the vortex to drive the TS 
wave but, as already mentioned, this interaction is of lesser importance than the 
mean flow modification due to the presence of the vortex. To simplify the analysis 
below we will hence indicate the sizes of these ‘mixed’ modes in the various layers 
but shall not need to give the full details of these modes. 

Firstly, we again concentrate on the main boundary layer where z = and 
we modify expansions (2.7 a - d) to consider 

« = r„u(0 + S ((Cf 10 + ...)£ + c.c.] + S 2 [r n e- 12 (t/ 30 + 

+ A [(£oi + • • -)F + 0{&€~ s E ±l F) + < 2 e- 12 (t> 21 F + . . .) + c.c] + . . . , (3 ' 2o) 

v = r nt}(£) + 6 [(Fio + ...)£ + C.C.} + S 2 [r„e“ 12 (F 3 o + ...)+...] + . . . 

+ A [(V„, +...)F + 0(Sr*E±F) +6h-' 2 (V 21 F +...) + c.c] + . . . , (3 ' 26) 
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(3.2c) 


w = c 12 tZ)(0 + 6 [(W 10 + ...)E + c.c.} + 6 2 [(W 30 + ...)+•••] + ••• + A 
[e 3 {W 01 + ...)F + O{6e- 3 E ±1 F)+6 2 e~ 9 (w 21 F + ...)+ c.c] + ..., 

p = € 24 p(£) + 5 [(P 10 .. .)£■ + c.c.] + 6 2 [(Pao + ...)+•■■] + ••• + A 

[e 3 (Poi + • • -)F + 0(8e~ z E ±1 F) + 6 2 e~ 9 (P 2l F + ...) + c.c] + . . . . 

Here the vortex terms (Ui o, . . . , P\o) and mean flow corrections (t/ 30 , • • • , P 30 ) 
satisfy the equations previously given in (2.8), (2.11). Inserting (3.2) in the govern- 
ing equations (2.5) leads to the usual TS solutions 

Uoi = Aou , Voi = Aov , Woi = -Uob, Poi = const., (3.3) 

r n 

where Ao denotes an unknown displacement and we have defined Ujs = r n 6tjU+PjV. 
We also find that 


C /21 — AoU 30 + Bqu , V 21 — ^ 0^30 + &ov > 

W 21 = — — — (Uoob + r n d 0 C /30 + PoVso'] ~ ~ — —Uq b, P 21 = const., (3.4) 

r n V ) r n 

where Bq is another constant and Uoob = r n^oo^ + Poo^- 

For the present we assume that these solutions are valid for £ > £; i.e. above 
the critical layer. An explicit analysis of the critical layer is required to determine 
how these solutions need to be modified for £ < £. 

Moving into the upper deck where z — e 9 y and y=0(l) say, the vortex decays 
exponentially and the expansions are 

« = A [e 3 (fi 0 + • ■ ■) F + S 2 e~ 9 (u 1 + ...)F+ c.c.] , (3.5a) 

v = — r„ + . . . + A [e 3 (6 0 + . . .) F + 6 2 e~ 9 (d x + . . .) F + c.c.] , (3.56) 

w =e l2 w(oc) + 6 2 W 3 o(oo) + . . . + (3 5) 

A [e 3 (& 0 + ...) F + S 2 e~ 9 (tZ>i + ...) F + c.c.] , 

p =e 24 p(oo) + A [e 3 (p 0 + . . .) F + <5 2 f~ 9 (p x + . . .) F + c.c.] . (3.5d) 

Feeding these expansions into (2.5) yields a straightforward upper-deck problem, 
see Smith (19796) for example. We demand that all disturbance quantities decay 
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as y — * oo and obtain matching conditions between the main deck and upper deck 
solutions. In particular, we find that 


•/« £. , Poa , . du>0 _ . n i. - - Po .p, » dp 0 

*(a 0 uo H t>o) + -x— = 0, -i(3 0 u 0 = -ia 0 p 0 , v o = — , -t^tuo = — rr' 

r n oy r n dy 


We need solutions of this system which match with (3.3) as £ — ► oo, and this 
requires 

f>2 

(3.6a) 


p _ Pi i 

“01 ~ A)) 

Tnl 0 


where 7o — &o + 7 ?-. Considering the next order system, matching with (3.4) is 
achieved if 


P *21 


A 

r n lo 


■B 0 = 


ZPoQlo ~ Po ^«o«oo + %f 00 ) 


Polo 


01 - 


(3.66) 


We need next to consider the critical layer where z = c 12 £ + e 16 z, with 2 = 0 ( 1 ). 
The main deck solutions (3.3) & (3.4) mean that inside the critical layer our fluid 
quantities develop as 


u =r n («£ + u'ge 4 z + ...) + 6 [(e 4 U 10 + ...)E + c.c .] + 6 2 [r n e 12 £7 30 + ...]+... 
+ A [(#oi +...)F + 0(6 2 e~ 8 E ±1 F) + 6 2 (<r 16 # 21 F + ...) + c.c.] , 

(3.7 a) 

v = r n(®f + v^e 4 z + ...) + 6 [(e~ 4 F 10 + ...)E + c.c. ] + 6 2 [r n t~ l2 V Z0 + ...]+... 
+ A [(#01 + . ..)F + 0{6 2 e~ 8 E ±1 F) + 6 2 {e~ 16 t 21 F + ...)+ c.c.] , 

(3.76) 

w + . . .) + 8 [(t^io + . . >)E + c.c.] + <5 2 [W 30 + . . .] + . . . 

+ A [e 3 (^oi +...)F+ 0(6 2 e~ 3 E ±l F) + 6 2 {e~ 9 & 21 F + ...) + c.c.] , 

(3.7c) 


p — c 24 p(£) + . . . + 6 [(P 10 + . • • )E + c.c.j 4- 6 2 [P 30 + . . .] + • . . 


+ A 


e 3 (#oi + ...)F + 0{6 2 e~ 3 E ±l F) + 6 2 {e~ 9 P 21 F +...) + c.c. . 


(3.7d) 
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The vortex terms U l0 , V 10 , U 30 ,. • ■ are as in expansions (2.15a - d) and fol- 
low from the critical layer analysis of the previous section. On using (3.7) in the 
governing equations (2.5) we find that the TS quantities are given by 

U 0l =A 0 n v Voi = A<4 # 01 = -^(£>ob) { ., (3.8a) 


U 21 —A(y 


dU; 


30 


dz ’ 


t> _ i ^ 30 
V 21 — Aq ~ — , 


W 2i = ~~ *Aq 


a 0 U 30 + —^30 

r n 


dz 
+ const.. 


(3.86) 


P 21 = Pil- 


Matching these solutions with the main deck quantities (3.3), (3.4) (which we 
supposed to be valid for £ > £) shows that, for 5 — ► -oo, then (3.8) matches with 
(3.3) & (3.4) if these latter solutions are the main deck solutions for £ < £ as well. 
Hence, as far as the TS wave is concerned, the critical layer at £ = £ is a purely 
passive affair and (3.3) & (3.4) continue to be the valid main deck solutions below 
the critical layer. To complete the analysis for the TS wave we need to consider 
the thin viscous wall layer which is of thickness 0(c 15 ). We notice that this TS 
lower deck is asymptotically much thicker than the 0(c 16 ) sized wall layer for the 
vortex. If inside the lower deck we have z = e 15 Z with Z - O(l), then the vortex 
quantities and the mean flow correction generated by the vortex are easily found by 
taking appropriate Taylor expansions about £ = 0 for the vortex terms satisfying 
(2.8) and (2.11). In particular, in this wall layer the fluid quantities take the forms 


u =r n u 0 e 3 Z + . . . + 6 [(^k> + • • •)& + c - c -] + ^ [ r * e 6 rn Z + • • •] + • ■ 
+ A [ (t4 + ...JF + 0{6e~ 3 E ±l F) + 6 2 {e~ 12 ul 1 + -~)F + e.c] + 

v =r n v 0 e 3 Z + ... + 6 (V^> + ...)!? + c.c.j +6 

+ A + . . .)F + 0(6e~ 3 E ±1 F) + S 2 (e 12 + . • • )F + c.cj + 

r 9g z 2 i 

e 3 (1^J 0 + ...)F + c.c.] + 6 2 e G w " m - — + -..J +••• 


(3.9a) 


9 v m Z + . . .j + . . . 


n 6 


W =W n £ 15 Z -t- . . . + S 


+ 


A [e 6 (W^ + . ..)F + 0(6E ±l F) + 6 2 (e~ G wl + ...)F + c.c] + . . . , 


(3.96) 

(3.9c) 
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V = £ 24 p(0) + ... + « [«*(*& + ••■)«+ C.C.] + 6 2 |p m (0) + . . .) + . . . 

+ ^ [‘Vo 1 . + • • -)F + 0(iE ±l F) + S 2 ((~ e pl +...)F + c.c] + . . . . 

(3.9d) 

Here u' 0 ,v' 0 ,w' 0 denote the values of ( fi ‘(fl) €=0 . («'(£)) , (©'(£)) ^ where 

u,v,w satisfy (2.3), and similarly u' m ,v' mi w„ denote the values of the derivatives 
of the mean flow terms 030, ^30 and W 30 (defined by (2.11), (2.13), (2.14) & (2.21)) 
evaluated at £ = 0. 

Substituting these expressions into (2.5) and comparing leading order terms in 
the TS disturbance quantities leads to the usual lower deck problem 


• / « , A) T ,t, dW^ 

1 { oc oUq 1 + — V^) + _|i = 0, 


i{X 0 Z - Cl 0 )uh + r n u' 0 W^ - -»a 0 pji + ~ JiV ~ , 

oZ* 

i(XoZ - ho)V^ + r n v' 0 wh = + 92 


dZ 2 ’ 


4 = >01, 


(3.10a) 

(3.106) 

(3.10c) 

(3.10d) 


with the necessary no-slip conditions 0* = = 0 at Z = 0 and the 

requirement that the solutions should match with the main deck solutions (3.3) as 

Z + 00 * Here we ^ ave denoted A 0 = cc 0 r n u' 0 + j3 0 v' 0 and we solve (3.10) in the 
usual way. The unknowns £/qi > ^01 are eliminated between (3.10a - c) to yield an 
equation for W which may be solved analytically by making the substitutions 


A* = tAo, n = A j ( Z - 2^- 

A 0 


no = 


a Jn 0 

a - 
^0 


(3.11) 


With these definitions it is found that the no-slip conditions at Z = 0 together with 
the matching with (3.3) can only be achieved if 


rnllPoi Ai(t)dtj = — t A 0 A f AqAi ( 770 ) , 


(3.12a) 


where At is the Airy function (see Abramowitz & Stegun 1964). Combining (3.12a) 
with (3.6a) yields the leading order eigenrelation, 

fn Ai W dt 


6 

U 3 


At (no) 


1000 


(3.126) 


16 



(3.13) 


It is well known that (3.12 b) has solutions for real &o,$o and Ho with 


= K! = 2.297, 


\ 3 
A o 


0 g% 


K 2 = 0.4355, 


see 


Smith (1979a, 6). Since our TS waves are three-dimensional, (3.13) admits an 
infinity of solutions. It is convenient to eliminate the radial distance r n by writing 


j ^ 3 A 1^ a __ 2. A 

fld = rn* n 0l a<i = r^do, & = r» T # o. Id = and A d = r n 4 A 0 . 


Then if (J. = 


@4- we obtain the solutions 

C«d 


4 1.001 (u 0 + /it/Q 

= 

Here n measures the angle <f> which the TS wave makes with the radial direction 
on the disc, since <f> = tan -1 ft. The solutions & d ,p d and fl d are sketched in Figure 
(2). The solutions (3.14) are, strictly speaking, valid only for fx < fx c (= 0.8284) or, 
equivalently, for X d > 0. The analysis outlined above has implicitly taken \ d > 0 
although the modifications required for X d < 0 are quite straightforward and stem 
from the consideration of the appropriate branches of the many-valued functions 
which arise. 

We may derive equations satisfied by u\ i , by equating higher 

order terms when (3.9) is substituted in (2.5). We obtain an equation for by 
eliminating u\ x & V 2 \ using a suitable combination of the r- and 0- momentum 
equations and the continuity equation, and we find that 


Pd - H&d, 


Q d = 2.297 of (uq + /i©o) • (3.14) 


= [n + (4 _ 
dr ] 4 dr] 2 A 2 Ai (r?o) 

a $ n 00 


rt 


(3.15) 


where A 0 o = ct 00 r n u 0 + fioov 0 , A m — a 0 r n u m + fiov m , $ — (% 00 '+X™)n 0 ’ an< * t ^ ie 

substitutions (3.11) have been made. The solution of (3.15) which satisfies the no- 
slip conditions at Z = 0 and which matches with solutions (3.4) in the main deck 
provides a second relation between P 21 and Bq. Combining with the earlier result 
(3.66) and using (3.6a), (3.12a) provides the eigenrelation for the correction terms 


17 


“005^00 and n 00 ;- 


t'A 0 


7oA 


ocq &oo + 


ftofioo 


t(Aoo + A m) ^ 2 
“7o 


+ 


^At(i7o) f_2K 1 X[ 
At (tj q) l 


+ 


*8 

Aq ^00 


2A 0 


g-r(Aoo + A m ) -f- 

0o 


^ (Aoo + A m ) 
2iXolo0oo 


+ 


3K 2 n 0 

2Kf\ 0 


Xq Kifloo 


3K 2 {Iq (Aqo + A m ) 


(3.16) 


01 


A A 

This equation admits real solutions for aoo,/?oo and ftoo only if the coefficient 

J 1 . A 2 . A. 

of Ai{r} 0 )/Ai (t? 0 ) vanishes. Redefining d 0 o = r n * ao 0d and 0 O o = r£0 ood implies 
that 

dood = M\0oo d + M 2 , floo = fn |i2i/?ood + j (3.17a, b) 

where 


Mi = 


and 


( 


2^ 


djjl 



3/? 2 


/ 5^u' n 

V W 2 d 





i2, = 


2ifi(M 1 UQ + Vq) 

3 { Q:dtt 0 ”1" 0d.v o) 5 


Mo = 


5 Id {adUrn + 0dVm ) 

3 0 2 d ( B \ Q 'o _ Ai&i 
V 30* w! 


)’ 


R 2 = 


2-K - ! + M 2 Uq) 

3(dd«o + ^d®o)^ 


(3.17c, d) 


(3.17 e,/) 


The dependences of M\,M 2 ,R\ & i?2 upon the parameter fi are illustrated in 

A 

Figure (3). We see that depending on the chosen value of 0ood, the neutral values 
for doo<£ an d flood can be of either sign. The results (3.17) indicate how the neutral 
curve for the TS wave is affected by the presence of the vortex and in the following 
section we demonstrate that the effect of the vortex is largely determined by the 
inclination of the TS wave (tan -1 fu) to the radial direction. In addition, we now 
consider the more realistic problem of allowing for a larger TS wave and indeed it 
is taken to be sufficiently large so as to be nonlinear itself. 
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4. The interaction of a linear vortex with a weakly nonlinear TS wave 

Our aim is to increase the size of the TS wave until the mean flow correction 
produced by the self-interaction of the TS wave is as large as the mean flow cor- 
rection produced by the vortex. This will ensure that at this stage the TS wave 
becomes nonlinear and this approach can be used to determine the stability of the 
TS wave to the vortex. We still assume that the vortex is small so that 6 C 1 and 
then analysis of the system described in Section 3 suggests that the crucial size of 
the TS wave when nonlinear effects become important is when 

A = <5<T 3 , (4- 1 ) 


where A is the unsealed amplitude of the TS wave defined in (3.2). It is straightfor- 
ward to check that at this TS size the analysis presented in Section 2 to determine 
the structure of the stationary vortex is not affected by interactions of the TS wave 
with either itself or the vortex. Consequently, to the orders required here the results 
of Section 2 & 3 for the vortex quantities continue to hold. 

As is usual in the weakly nonlinear problem for a TS wave the nonlinearity 
manifests itself in the lower deck, whereas the main and upper deck solutions remain 
essentially linear in character. We slightly modify the method used in the previous 
section by obtaining an evolution equation for the TS amplitude in the vicinity of 
a chosen point rather than obtaining an equation for the neutral curve and we base 
our approach on that given by Hall & Smith (1982). If we consider TS waves close 
to the point (r n ,0 n ) we define the O(l) coordinate r x by 

r = r n + 5 2 e- 12 r 1 , (4-2) 


and allow the TS wave to evolve on the r x lengthscale. Consequently, we redefine 
F given by (3.1) by 


F = exp 




0(3{e) - n(c)t 



(4.3a) 


where, 


a =oco, & = 0o + •••)+ 12 (/?oo + ■••)» 

n = e 3 (n 0 + . . .) + s 2 e 9 {h 0 o + ...)• 


(4.36) 
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Here the vortex terms u} 0 ,...,P* 0 again follow from suitable Taylor expan- 
sions of the solutions to (2.8) around £ = 0. A significant difference between the 
expansions (4.5) and the corresponding (3.9) for the linear TS problem is that now 
the mean flow correction induced by the vortex in the lower deck is of the same 
size as that produced by the interaction of the TS wave with itself. Hence the 
terms U mL ,. . . , P mL are no longer the simple Taylor expansions of the mean flows 
(U 30 ,...,P 30 ) generated by the vortex interacting with itself in the main bound- 
ary layer and which satisfy (2.11). Since outside the boundary layer the TS wave 
self-interacts to produce mean flow terms of size only 0(6 2 e -9 ), (i.e. the same size 
correction as in the lower deck) we notice that equations (2.11) are still relevant, 
together with boundary conditions (2.13), (2.14) and jump conditions (2.21). Hence 
the leading order mean-flow terms across the whole boundary layer are unaffected 
by the presence of the weakly nonlinear TS wave. 

Substituting (4.5) into equations (2.5) and examining the TS quantities yields 
a system which is very similar in character to that discussed by Smith (19796). At 
leading orders we find that the linear lower deck problem (3.10) applies so that 
matching with the main deck solutions yields the linear eigenproblem (3.126) as 
before which has solution given by (3.14). 

The mean flow terms U mL and V mL satisfy equations of the forms 


r n 


d^UrnL 

dZ 2 



-4(4r) + wlir^f- 


+ W C 


t awfi) - 


01 


dZ 


(4.6a) 


d 2 V, 


dZ 2 


ta ° (i U 0 l)* V 0 l Uoi( v oi)*) + ( W oiY~^r + W 01 ~ Q 2 ^ • (4.66) 

To match with the mean flow produced by the self-interaction of the vortex in 
the main deck requires 


U, 


mL 


VmL * V m Z y 


as 


00 , 


where u m and v m are defined below (3.9). Substituting the solutions of (3.10) 
together with definitions (3.11) reveals that U mL and V mL are given by 


UmL = 




t 



(4.7a) 
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V mL = -4t- ( f a M<“) Hoi 2 + -%r(» - -Jo), (4.76) 

* ^ Tn A j VAo / *3 Aj 


where 


Gi(«) =r» 


C 2 («) = i -7 


2 rt 

/3d(/?ci^Q ~ &dV 'o)^(t) jrry r / Ai(r)cir 

f Ai(r)dr 
•'no 


| A / 

ctd{ocdV 0 - &fU 0 K(*) - 


r = 


Ai'(r t 70 ) 

fno (ti Mrjdr) dt, 
Ai'(rio) 


T* + c.c., 


r + c.c., 


and where £(t) satisfies — t£ = 1 with £(r?o) = 0, £ — ♦ 0 as t — ♦ oo, so that 

r* (/£ -A*( r ) rfr ) 

Uo Mh*i)) 


The second harmonic terms in (4.5) satisfy the equations 

2>(So Ul + ^4) + = 0, 


(4.8a) 


2i(A 0 Z - fio)t4 + r»aiw , J t , = -2iao^, + oj, - W„ t , ^ + ^f 1 , (4.86) 

2i(A 0 Z - f) 0 )4 + = - ^4, + v oi - w m (4-8c) 


Pj X = const. 


( 4.8<f) 

We may obtain an equation for by eliminating £7^ and from (4.8a, 6, c). 
We require the boundary conditions of u \ x , and vanishing on Z = 0 and 
matching conditions as j? — * oo are obtained by considering the main and upper 
deck equations to be satisfied by the harmonic terms. This is all straightforward 
and follows from the work of Smith (19796) so we merely state the solution for W \ x 
here. We obtain 


d 2 W. 

1 hj 1 


t 


21 _ 


2iAo(.-4o) 2 

r l (f™Mt) dt ) 


rTJ 1 

F* + At (rj) I Ai(t)dt + B^Ai(2^rj), 
Jrio J 
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(4.76) 


VmL ~ Trh 3 ( / g 2(0^) l^o| 2 + —TZ I (v - Tjo), 

\Jt)o J t'ar^AJ 


where 


Gi{t) = i~* 


G 2 {t) = i~* 


2 /*t 

PdiPdU o - du<>o)<£(*) - > / At(r)dT 

A * fao) y no 

# a / J 2 rt 

a<*(“<*®o - AlfioK(t) - ~ A V ?? d x / At'(r)d7 

(»7o) 

fno (C o ^*( T ) rfr ) 


r* + c.c., 


r* + c.c., 


r = 


Ai' ( t] 0 ) 

and where £(t) satisfies - it - 1 with £(*7 0 ) = 0,.C 


0 as t — ► oo, so that 


f 

ho {Ax{t i)) 2 


The second harmonic terms in (4.5) satisfy the equations 

A 

±v, T ,) + ^ 

r n 


2,-(&o4 + ^4) + ^ = 0, 

2i(X 0 Z - h o)4 + r „s>t = _ 2ifiB pt + 4 *4. _ art + 4^|L, (4 . 

t 


2i(X 0 Z - fl 0 )V 2 \ + r^K = -22 J»£ + 4 ^ 


t - _ 2l ^° Dt , T/t 5W 01 H /f 5V 01 . " '21 /, o 

^01-^-+ > ( 4 - 8c ) 


dZ dZ 2 

rt 

• + 


(4.8a) 

86) 


dZ 


dZ 2 


P^i = const. 


(4.8d) 


We may obtain an equation for W £ by eliminating u \ , and v£ from (4.8a, 6, c). 
We require the boundary conditions of 17^,1^ and W 2 \ vanishing onZ = 0 and 
matching conditions as Z > oo are obtained by considering the main and upper 
deck equations to be satisfied by the harmonic terms. This is all straightforward 
and follows from the work of Smith (19796) so we merely state the solution for wl 
here. We obtain 

d 2 W^ 


dr)' 


— 2lA °( yto ) — 4 + Ai’(ri) [ Ai(t)dl + B t Ai(2*‘f ) ), 

r »(/“ A, •(»)*) L J "° J 
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where 


and 


e* ) [/* At (92) £(42)^92] 

F* = Ai(rj) / 1 777 772 

Jfjo (-At(flx)) 

F = -2 - * ]^2Ai" {T]) Ai(rj) + Ai {r)o)Ai (»?)] , 


dqi. 


i | = 2 *r), 


and the constant 5| is chosen so that dW 21 /dr) - 0 at tj r? 0 and 



A 2 

A o 

a 3 ^ 

dr] 

A* 1 oPo& $ ■ 

dr ] 3 


as t] 


00. 


(4.9) 


ri-Vo 


Finally, to recover the evolution equation for the amplitude parameter A 0 for 
the TS solutions we inspect the governing system for the terms U 31 ,V 31 ,W 31 ,P 3V 


This set of equations take the form 


A 


■( - rrt , ^ , dW Sl _ p 

t(otoU 3i + ^ ’'31) ^ Qg 


i(\oZ — ^0)^31 + r n ti 0^ r 31 * **0-^31 dZ 3 F 2 ’ 

A + , t id 0 t d 2 V^ 

i{XoZ — + r n v 0 W^ x 4 


it a 2 14, _ 


pi _ “ ?-L = Ft, 

*31 aZ 2 

pJi = const., 

where F X ,F 2 and F 3 are combinations of lower order TS and mean flow terms. We 
can use the standard method to obtain a governing equation for W 3l which assumes 

the form 


d 3 W. 


•t 

31 


dr ] 3 


- V 


a < + wl = ^Bo + ^ 


drj 


r n A, 


r n A* 


2/3 00 4. 1 f i«n — + 

A> Tfo V dTl n 




A °- a 7 


„,t i’K 1 ,)' latviao*?,)* i w t r £!l 4 

az } - - 5 -37 37^ 2' 01 1 3Z 2 


2 az az 


^i^L(a 0 r n U m L + PoV mL ) + *wJi -^{& 0 r n U mL + 
oZ c ^ 


/K A 

(Aoo + A m) 
tA* 


(fo - Vo) + *Vo) - Wji 


(4.11a) 
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where we have used definitions (3.11) and those given below (3.15). The relevant 
boundary conditions are that 


^1 = 0 = 


dW. 


31 


dr] 


on rj = rjo. 


and 


as rj 


oo. 


dW. 


31 




dr) 


7 /\ * 

r n*-* * 


(Aoo + A m) 


iB 0 Xo 


r n A: 


(4.116) 


The system (4.11) only admits a satisfactory solution if a certain solvability 
criterion is met, and it is this criterion which leads to the desired evolution equation. 
To obtain this condition, we need to consider the adjoint problem to (4.11), and 
using the results of Ince (1956), Hall & Smith (1982), it is found that the adjoint 
function for this problem is 


S(r!) = At{ rj) - (4.12) 

where the function £ is as given in the definition for and G 2 in (4.7). If we write 
the right hand side of (4.11a) as R*(r}), multiply both sides of (4.11a) by S(rj), 
integrate by parts and use (4.116), we obtain the solvability condition 

iPlloAi{rio) ( . d A ,\ , f°° 

r n A? ^ ir nU 0 ~ + ^oov 0 jAo = J SR*dr}. (4.13) 

We now use the definitions of i?* together with the solutions for U mL ,V mL and 

W m £ (given in (4. 7), (4. 9)) and that for wjj to obtain an evolution equation for Ao. 
This equation takes the general form 


Ti&oo + r 2 n 0 o + 73 ] ^4o + '-^-Aq\Ao\ 2 , (4.14) 

where T U ...,T 4 are complex functions of a d , (3 d , i d and Vo . All of the integrations 
necessary in the determination of Ti,...,T 4 were evaluated numerically. Initial 
conditions at a suitably large value of 77, say ij <», were found asymptotically and 
then the defining equations for the various functions integrated using a fourth order 
Runge-Kutta scheme. The integrations were all performed using the trapezium rule 
and the results were checked by varying the value of r}^ and the step length used. 
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These checks lead us to believe that the computed values of Ti - T 4 are correct to 
within 0 . 1 %. 

In the above equation, each of the coefficients T x - T A is a function of the 
orientation of the TS wave and so are functions of the parameter fi. Further, T u 
T 2 & T 4 are all independent of vortex quantities; i.e. are functions only of the TS 
wavenumbers a d J d and r? 0 . Since our prime concern at the outset was to determine 
the effect of the vortex upon the stability of the TS wave, the precise values of Ti, 
T 2 and T 4 are not of immediate interest, save to remark that we found that for all 
admissible values of (< Me ~ 0.8284), Re(T 4 ) < 0, see Figure (4). 

If we write (4.14) in the modified form 


^- = b<A 0 + cJo\Ao\ 2 , ( 4 - 15 ) 

dr i 

then it is clear that on the basis of linear theory the TS wave is unstable if Re{K) > 
0. However, (4.15) admits an ‘equilibrium’ solution with 



5 


so for Re{c .) < 0 a non-zero, finite, steady solution is possible. It may be easily 
shown that this solution is stable and then the TS wave is said to be supercritical^ 

stable. 

The only term in (4.14) which contains vortex-induced quantities is the term 
T 3 which is given by the formula 


. / 5 7d u o 


2Ki. _• .s. f K x 

A 1 d AiM 

. V 3 fid 

*<dfc) 

3A| Uljl V 1.0°1 

) Ai'{rio) 


T z 


= — 1 


5t 

W d 


lard / us\ d 

"wi-a (-*- - 1 ) #4] (**. + 

3^1 ld \ 1.001 ) Ax (T 7 o) 


(4.16) 


Here, Ki is defined by (3.13), m by (3.11), the nondimensionalised TS wavenumbers 
are a d J d and u' m ,4 are the derivatives of the mean flow terms U 3 0 , V 30 induced 
by the vortex when evaluated at * = 0. We see immediately that, for a given n, 
Tz is proportional to the size of the mean flow correction due to the vortex and 
on examining equations (2.11), (2.13), (2.14) and the jump conditions (2.21) which 
determine this mean flow it is clear that u' m and v m are, in turn, proportional to 
the square of the amplitude of the crossflow vortex. 
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Hence the effect of the vortex on the stability of the TS wave can be measured 
directly through the coefficient T 3 . More particularly, we can see from (4.14) that 
if Re(T 3 ) < 0 then the TS wave is stabilised by the presence of the vortex, whereas 
the opposite conclusion may be drawn for Re(T 3 ) > 0. This function is illustrated 
in Figure (5). We notice that Re(T 3 ) may be of either sign and this is determined 
solely by the value of ji (or equivalently by the orientation of the TS wave). In 
particular, we find that for /x > 0.297 or -0.296 < n < 0 (corresponding to TS 
waves making angles between -16.5° and 0° or greater than about 16.5° with the 
outward radial direction) the effect of Re(T 3 ) is destabilising. Conversely, for TS 
waves making other angles with the radial direction, the opposite effect is observed. 

Crucially then, our analysis in this section suggests that on a weakly nonlinear 
basis the TS waves are supercritically stable as is the case in the absence of the 
vortex, but that the presence of the of vortex can be both of a stabilising or of 
a destabilising effect, with the change between these states occuring for TS waves 
inclined at roughly 16.5° with the outward radius on the disc. Some comments and 
brief conclusions concerning these results will be made in Section 6. 


5. The secondary instability of a fully nonlinear vortex to a TS wave 

In this section we attempt to generalise the work contained in Sections 2 & 

3. Thus far our concern has been only with the problem involving a linear vortex 
and now we consider the effects of nonlinearity of this instability mode. We can 
anticipate that the nonlinearity will affect the governing equations for the mean 
flow induced by the vortex and that this in turn will alter the evolution equation 
for the TS wave. 

For simplicity, we shall revert to considering the problem of the interaction 
involving an infinitesimally small TS wave (of size A as defined in (3.2)). In this 
case, the nonlinear vortex structure is unaffected by the presence of the TS wave, at 
least to the orders that we will be concerned with. We can use the results described 
by Gajjar (1989), who showed that when the vortex size reaches 0(R~s) (=0(e 5 * * 8 )), 
then the first effects of nonlinearity are encountered in the critical layer although 
the remainder of the flow structure remains linear. It is this sized vortex which we 
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consider here and we note that the mean flow correction generated by the vortex 
is now 0(e 4 ) which is larger than the vortex itself. This large mean flow correction 
is necessary due to the properties of the nonlinear critical layer, see Stewartson 
(1981) and Haberman (1972). Since the vortex induces an 0(e 4 ) correction and 
following the ideas presented in Section 3, we expect the TS wave to have spatial 
and temporal dependence as in the function F\, defined by 


F i 



d(r, e)dr + 60(e) 


n(e)i | 


(5.1a) 


where, 

a =& o, 0 — (0o + •••) + e *(Poo + •■ •)> , 

(5.16) 

Cl = £ 3 (Ho + ...)+ e 7 (noo + •••)■ 

Here, as in (4.3), 0 OO and f2 0 o denote the changes in the azimuthal wavenumber and 
frequency of the TS wave due to the presence of the vortex. We permit the scaled 
amplitude of the TS wave to evolve on an r\ lengthscale, where 


r = r n + e 4 n, 


(5.2) 


replaces the definition (4.2). The linear theory for the TS wave will allow us to 
derive an evolution equation of the type (4.14) with the nonlinear term on the right 
hand side of the equation absent. Indeed, following the work contained in Sections 
3 & 4 it is easily shown that the stability of the TS wave is governed by an analysis 
identical to that already performed and that the linear evolution equation of the 
TS disturbance will satisfy an equation of the form 


dA 

dr i 


Q 1 0OOd + <3 2 Ciood + Qsj -<4> 


(5.3) 


A A _ £ A ^ i. A 

where A is the scaled amplitude of the TS wave, 0ood — r n J /?oo> flood = r n 
and Q u Q 2 & Q 3 are functions of the orientation of the TS and certain mean flow 
correction parameters which follow from the analysis of the nonlinear critical layer. 
Further, as in the discussion following (4.14), it is only the function Q 3 which 
contains the mean flow velocities induced by the vortex and hence is the only term 
to give information concerning the effect of the vortex on the stability of the TS 
perturbation. Finally, Q$ is given by the same formula as determined T 3 save that 
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the mean flow derivative terms present on the right hand side of (4.16) should be 
replaced by the equivalent terms for the (larger) vortex-induced mean flow. 

To summarise the above, we conclude that the evolution equation for the TS 
wave follows immediately once we have determined the mean flow correction due to 
the nonlinear vortex. We consequently consider this problem, in a manner similar 
to that used in Section 2 and below we determine the governing equations for the 
mean flow. 


5.1 The nonlinear vortex 

We follow Gajjar (1989) and use scalings for the vortex appropriate to the 
nonlinear critical layer calculation to be addressed. In the inviscid zone, where 
z = e 12 £, the implied velocity and pressure fields take the forms 


**n^(£) "1" ^ r n Um 4” C {U 10 T U m j) "f" . . . , 

(5.4a) 

r nv(£) + e 4 r n v m + e 8 (Vi 0 + u m i) + . . . , 

(5.46) 

u, = e 8 Wio + e 12 tZ>(0 + e 16 iD m + . .. , 

(5.4c) 

P — e*Pio + • • • , 

(5.4tf) 


where u m , v m and u> m are the leading order mean flow terms induced by the vortex. 
We reiterate that all the mean flow terms in these expansions (terms with subscript 
m) are functions of £ and r j and that the mean flow corrections in the radial 
and azimuthal directions are larger than the fundamental vortex terms. We seek 
solutions for these latter functions in terms of the crossflow variable x, where x is 
defined so that (2.6) may be rewritten as 


E = exp(fx). (2.6) 

Substituting expansions (5.4) into (2.5) and comparing suitable terms yields 
that there is a solution for the vortex terms of the form 

{U lo ,Vio,W l0t Pio) = 

(5.5) 

M u oi{ r n, 0 cosx,t>oi(r n , 0 cosx, Woi(r n , f) sinx,p 0 i(r n , £) cos x), 
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where A is the scaled size of the vortex and (u 0 i, voi,«>oi>Poi) satisfies 


q 0 uoi + — Uoi =0 , -UqbUoi + r n u w 0 i = «oPoi, 

K ( 5 . 6 ) 

' 00 rr - ' 

-UobVoI + fnV VJoi = P 01 , ^OB^Ol ~ ~Poi- 

r n 

We recall that U 0B = r n a 0 u + fi 0 v and note that tv 0 i satisfies the Rayleigh equation 
(2.9). 

Additionally, we find that the mean flow terms satisfy equations very similar 
to (2.11), and in particular 

(5.7a) 


d\ Vyn 

2u m H — = 0, 


dt 


d 2 u 


m _ du m 


— w — — 2 uu m — u u) m + 2(1 + v)v m — 0, 

at 


d 2 v m dv 

— w 


m 2uv r 


v' w m - 2(1 + v)u m = 0. 


(5.76) 


(5.7c) 


dt 2 d£ 

The wall layer structure for this nonlinear vortex is analogous to that in Section 
2 for the linear vortex and so we conclude that we need boundary conditions for 
the induced mean flow of the form 


u m = v m = u>m = 0 on £ = 0, u m ,Vm — ► 0 as £ * oo. (5.8) 

The most dramatic change between the work here and that given earlier is 
the treatment of the critical layer zone. For now we expect a nonlinear analysis 
to become appropriate in order to compute the jumps in the mean flow and its’ 
derivative across £ — £• As before, we suppose that the critical layer is described 


by £ = £ + e 4 z,z = 0(1), and then the analysis of system (5.6) as £ — 
that in the critical layer the relevant expansions are 

->• £ implies 

u = r n ti| + e 4 (u 0 + r n v,£zj + e 8 ^ui + -r n u^S 2 ^ + . . . , 

(5.9a) 

v = r n U| + e 4 (vo + + e 8 ^t>i + “ r nP| + . . . , 

(5.96) 

w = e 8 u>o + e 12 (d>i + + . . . , 

(5.9c) 

p = e 8 Po + c 12 Pi + • • • • 

(5.9d) 
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Using these expansions in (2.5) reveals that the first order unknowns u>o,po in 
(5.9) are given by 


» • A%p 0 o . 

wq — tn* sin x = — — sm x. 


B, 


01 


Po — Apoo COS x. 


(5.10) 


where p 0 o = Poi(£) as defined by (5.6) and 2?y* = r n ctju( k ){£) + PjV^(^). It is 
clear that these solutions match on to (5.5) and (5.6) outside the critical layer. 

If we define the unknowns uy* = r n aytifc + /?yv*, then we find upon substitution 
in the continuity equation (2.5a) and on using (5.10), that tioo = Co o (a constant). 
This, when matching with the inviscid layer solutions surrounding the critical layer, 
gives that 

lim 
«— * fH 

or, equivalently, that 


_ , ‘ 
"t" Vm 

= lim 

- , /?o _ 

a 0Um H V m 

r n 

€—*i- 

r n 


MOlfc = 


r n OL o r _ 


[«m(f)]£ , 


p o l “ mvs ' J £- 

i.e. the jump in t» m across the critical layer is proportional to the jump in u m across 
that layer. 

After straightforward analysis (see Gajjar (1989) for further details) it is found 
that t) 0 satisfies the equation 


(t) - . d i , AioPoo . dv 0 

(B 0l z + B l0 )— + -g_smx— = Ap 00 


Po 

r n Bqi 


d 2 v 0 , 

SmX+ ~dj*' ^• 11 ^ 


where B i 0 = B i 0 + C 0 o- 

We recall that Boi < 0 and this enables us to write (5.11) in a canonical form. 
If we define 


v 0 = K 0 (Y-V*), C*y = |^ + 5, Ko = ^^( — - r -^^) > 


>01 


BoiC't = K l , (> 0 ) c 2 = 


% \ r n 

2 . 


v 2.-.' 

7 C 

B oi 


2 _ ^7oPOO 


R 2 ’ 
n oi 


we can reduce (5.11) to the form 


BV* . dV * % d 2 V* 

y — h smi— - A c -XT?;- = 0. 


dx 


dY 


dY 2 


(5.12) 


(5.13) 
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We note that C* < 0 so z — * ±oo corresponds to Y — * Too and the boundary 
conditions to ensure a match with the inviscid zone solutions now take the forms 


COSZ 

K 0 Y 


(y 2 ) 


(5.14) 


as Y — * dcoo where (v m )^ denotes lim^ ,^(*ro) respectively. 

The system (5.13) & (5.14) is a classical one and the solution characteristics are 
well known, see Haberman (1972), Stewartson (1981). If we define V* = d 2 xjj/dY 2 
then we obtain the system 


. d 3 0 

dxdY 2 


d 3 0 , 3 4 0 

+ S1I1I— rr = A C ' 


' dY 3 


dY 4 ’ 


(5.15a) 
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■4 -Y 2 - + In 171 cos x + B${x) + ..., as Y — * ±oo. (5.156) 

2 Jfo 


We may apply the standard technique of integrating (5.15a) with respect to Y, 
applying the boundary conditions (5.156) and then integrating over a period in x. 
After appealing to periodicity and a further integration with respect to Y we find 
(see Haberman 1972) that, on using (5.12), 

(®m ~ ®m) = “ 2^'“ J Q ( B 0 ~ B 0 ) S1X1 xdx * ^ 5 * 16 ^ 

The numerical solution of (5.15) is well documented (Haberman (1972), Smith 

A 

& Bodonyi (1982)) and if we define the phase shift (-0) by 


1 f 2n 

-<j> = - / (Bq - Bo) sin xdx, 

7T Jo 


(5.17) 


then the dependence of 0 on the parameter A c is shown in Figure (6). (This figure 
was plotted using data kindly supplied by Dr. J. Gajjar.) It is known that as A c 
becomes very large 0 — * it which corresponds to the result for a linear vortex, 
whereas for small A c , 0 — ► 0. Using the transformations (5.10, 5.12) we find that 
(5.16) becomes 


|a m (e)]fc 


A)(/W| - 


^I'llBoi 


( 0 ) ( u >*) 2 > 


MO] |! = 


r n ao 


Po 


MOlft- 


(5.18a) 

(5.186) 
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These are the mean flow jump conditions which are the counterparts of (2.21). 
We can easily show that for A c » 1 (and hence $ — ► n), the size of u>* diminishes 
(from (5.10), (5.12)) and the jump criteria approach those of (2.21) for the linear 
vortex. 

We need to extend the above to consider the shift in the derivative of the mean 
flow across the critical layer and this necessitates consideration of the next order 
terms in expansions (5.9). 

To derive the necessary information we can proceed to obtain a differential 
equation to determine tq. The governing equations for these unknowns are given 
by Gajjar (1989) so we do not repeat them here. Instead, we merely state that we 
can reduce the problem for Vi to the canonical form 


d 3 * d 3 x a 4 * , K 0 Wf 

Y dlw + smz W* - X ‘W* =x(l ’ r)-- cT’ (5 ' 19a) 

with boundary conditions 


y3 +~- [ J B 0 i(v ml ) qF - BxottO*] Y 2 +0{Y In |y|), 

(5.196) 

as |F| oo. Here (t4) T , (”mi) T denote [4(0] , [®mi (01<— respec- 

tively, where these mean flow terms are defined in the expansions (5.4). Also x(z, Y ) 
is a very complicated expression involving ip (the solution of (5.15)). However, we 
find upon integrating (5.19a) with respect to Y and integrating over a period in x, 
that the function x{ x >Y) plays no part in determining the jump in v' m across the 
critical layer and, on reapplying transformations (5.12), we obtain the result 

[®rn(0]!_ = • (5.20a) 

Similarly, we also find that 

[®m(0] f _ = ®f[6m(0]ft • (5.206) 

Consequently, we now have a complete determination of the mean flow induced 
by the nonlinear vortex; namely the defining equations (5.7), boundary conditions 
(5.8) and jump conditions (5.18) & (5.20). This system is almost identical to that 
studied in Section 2 and in practice (u m >t i m,u) m ) is merely a multiple of the so- 
lutions (U 30 ,V 30 ,W 30 ) found earlier. This multiplication factor arises due to the 




i 

6 


1 + r n B 0 iC?{v' m )* 


32 



difference in the jump conditions (5.18a) and (2.21a). This observation leads us to 
some immediate conclusions concerning the effect of the vortex on the stability of 
an infinitesimal TS wave. For a given orientation of the TS wave, the function Qz 
appearing in (5.3) is a real multiple of T 3 in the evolution equation (4.14). Conse- 
quently, the nonlinear vortex has a stabilising influence on the TS wave if the latter 
mode makes an angle <f> with the outward radial direction if 0 < <j> < 16.5 or if 
(j> < -16.5°, and is destabilising otherwise. 

To conclude this investigation of the weakly nonlinear vortex problem we con- 
sider the size of the mean flow jump across the critical layer. We know from the 
forms of (5.18a) and (4.16) that the stabilising or destabilising effect on the TS 
wave is proportional to the jump across the critical layer. From the transforma- 
tions (5.12), solutions (5.10) and (5.18a), it is seen that the jump in u m across the 
critical layer is proportional to <f>A which is a function of A c . Figure (7) illustrates 
this dependence, and two asymptotic cases naturally arise. Firstly, as A ► 0 we 
know that A c — ► oo and 4> — + n so the jump across the critical layer becomes 
small. This corresponds to returning to the linear theory of Section 2. Addition- 
ally, for small A c (which corresponds to a large vortex), we have from (5.12) that 

2 

the scaled amplitude of the vortex A ~ 0(A C *)■ It is also known that in this limit 

v A c C (1) , where = 1.379 (Smith & Bodonyi 1982). Hence the mean flow 

jump is 0(A7® ) and so becomes large. Following the work of Gajjar (1989) it can 
be shown that when the scaled amplitude A ~ 0(e - ~), the nature of the critical 
layer changes to a structure similar to that given by Bodonyi et al (1983) for a 
strongly nonlinear critical layer and thus a modified analysis of the flow is required. 

We have thus shown in this section that the effect of a nonlinear vortex on 
the stability of the TS wave is very similar to that of the linear vortex considered 
earlier with the main difference arising from the changed mean flow jump across 
the critical layer. As previously, the nature of the influence of the vortex is largely 
determined by the orientation of the TS wave. Of interest would be the extension 
of our work to study the case of a strongly nonlinear vortex; the flow structure for 
this problem following from the discussion of the limit A c ♦ 0 given above. 
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6. Conclusions 


In this paper we have attempted to provide a rational asymptotic analysis of 
the problem of interaction between a stationary crossflow vortex in the flow induced 
by a rotating disc and a classical lower branch TS wave. The interaction has been 
studied in the context of this particular basic flow because this flow is susceptible 
to instabilities which occur in the boundary layer of a swept wing and this is of 
relevance to the development of Laminar Flow Control wings. 

At first, the asymptotic structure of a linear crossflow vortex was obtained and 
we have investigated the effect of this vortex on both linear and on weakly nonlinear 
TS waves. Of particular interest in our work has been the problem of the behaviour 
of the TS wave in the neighbourhood of the critical layer of the vortex which is 
situated inside the main part of the boundary layer of the flow. As far as the TS 
wave is concerned, the critical layer plays a passive role and the leading order TS 
solutions pass through the critical layer region unaffected. In addition, we have 
demonstrated that the effect of the vortex on the stability of the travelling TS wave 
is felt entirely through the mean flow generated by the presence of the vortex. On 
analysing this phenomenon, we have found that depending on the orientation of 
the TS wave, the vortex can have either a stabilising or a destabilising effect- in 
particular the vortex stabilises the TS wave if the latter makes an angle between 0° 
and 16.5° or less than —16.5° with the outward radial direction. This work has an 
obvious practical implication, namely that it is possible that when these two types 
of instability are present in a three-dimensional flow that the crossflow vortex can 
destabilise the TS wave. 

The above interaction structure has also been extended to the study of a non- 
linear vortex. The findings described in the last paragraph are largely unaltered 
by this change, at least as far as the effect of the vortex on the stability of an 
infinitesimal TS wave is concerned. The analysis for the linear vortex was easily 
adapted to study the implications for the stability of a weakly nonlinear TS wave in 
the presence of a linear vortex, although the work for the case of a nonlinear vortex 
is not so readily extendable for the problem with this larger TS amplitude. This 
extension is currently being considered by the authors. 

Further, the study of Section 5 and that of Gajjar (1989) points to an even 
larger crossflow vortex which may be described by asymptotic means. This again 
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provides scope for future work. 

Finally, we should remark that here we have considered the problem of inter- 
action between a pair of specific crossflow and TS instabilities. In realistic flows it 
is likely that other modes could be present and it is desirable to be able to classify 
the relative importance of other possible interaction mechanisms. In particular, 
we are looking at the problem of interaction between TS waves and the stationary 
vortex mode described by Hall (1986) which is characterised by having zero wall 
shear stress for the effective crossflow velocity profile. In addition, we note that our 
approach developed here can only deal with interactions at asymptotically large 
Reynolds numbers and the importance of the crossflow vortex on the stability of 
the TS wave at lower Reynolds numbers can only be resolved by pursuing extensive 
numerical calculations. 
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Figure Captions 


Fig (l). The dependence of the mean flow terms U 3 o,V 3 o and W30 (defined by 
(2.7) and satisfying (2.11), (2.13), (2.14), (2.21)) upon the boundary layer coordi- 
nate £. 

Fig (2). The neutral non-dimensional wavenumbers aj, Pd and frequency f 2d for 
the TS wave defined by (3.1) expressed as functions of the waveangle <f> = tan 1 /x. 

Fig (3). The functions (i) M U M 2 and (ii) Ri,R 2 which determine the correc- 
tions to the neutral wavenumbers and frequency of the TS wave due to the presence 
of the vortex. 

Fig (4). The real part of the coefficient T 4 of the nonlinear term in the TS 
evolution equation (4.14) expressed as a function of /x. 

Fig (5). The real part of the coefficient T 3 in the evolution equation (4.14) 
which determines the effect of the crossflow vortex on the stability of the TS wave. 


Fig (6). The function $ = 4>{ A c ) given by (5.15) and (5.17). 

Fig (7)- Dependence of the vortex quantity 4>A 2 upon the parameter A c . Here 
4> is defined by (5.15), (5.17) and A is the scaled vortex amplitude size. 
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